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Chapter  I 
INTEODUCTION 

In  the  late  1920's  a  nnmber  of  experiments  were  condneted  to 
investigate  the  properties  of  high  intensity,  ultrasonic  vibrations  in 
liquids.  The  results  led  to  similar  experiments  in  gases  where  it  was 
noted  that  particula  '.es  suspended  in  the  gas  tended  to  cluster  and  form 
aggregates.  Thus  began  research  in  the  field  of  acoustic  (or  sonic) 
agglomeration  (or  coagulation).  As  described  by  Mednikov  (1965)  perhaps 
the  most  notable  of  the  early  experiments  in  this  field  were  those  of 
Brandt,  Freund,  and  Hiedeman  using  tobacco  smoke  suspended  in  air. 
Among  other  things  they  showed  the  strong  frequency  dependence  of  the 
phenomenon.  In  particular  they  discovered  that  agglomeration  occurred 
more  readily  when  sound  frequencies  were  reduced  to  within  the  audible 
range . 

Through  the  1930's  the  interest  in  exploring  the  fundamental  nature 
of  the  subject  waned  while  the  search  for  practical  applications 
continued.  According  to  Boucher  (1960)  it  was  in  1931  that  interest  in 
the  meteorological  application  of  the  process  first  occured  when  A.  Amy 
took  out  a  U.  S.  patent  for  a  device  to  clear  fogs  around  airi'Orts  by 
agglomeration  of  droplets  to  precipitable  size.  The  device  war  to 
consist  of  a  set  of  whistles  arranged  in  a  gimbaled  frame  which  was  to 
be  attached  to  an  aircraft.  When  the  aircraft  flew  over  or  through  the 
fog  the  device  would  emit  a  sonic  beam  which  could  be  swept  back  and 
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fc*th  like  a  searchlight  and  saturate  the  clond  vith  acoustic  energy. 
It  is  not  mentioned  whether  the  device  was  actually  ever  assembled  or 
tested . 

Significant  advances  in  the  application  of  sonic  agglomeration  cjme 
during  the  World  War  II  years  when  efficient,  high  energy  output  sound 
sources  became  available  in  the  form  of  sirens  manufactured  by  C.S. 
firms.  Since  then  there  hns  been  continuous  progress  in  the  area  of 
industrial  application.  Uses  include  the  removal  of  particulate  matter 
from  the  exhaust  of  industrial  plants.  It  has  been  proposed  (Shaw  and 
Rajendraa,  1979)  that  acoustic  agglomeration  could  be  used  to  control 
the  radioactive  sodium  particles  that  would  result  from  a  fire  in  a 
Liquid  Metal  Fast  Breeder  Reactor  plant. 

There  have  also  been  some  attempts  at  applications  of  acoustic 
agglomeration  to  the  atmosphere.  Most  of  the  work  was  done  from  World 
War  II  through  the  mid  1950's  and  was  directed  toward  clearing  fog  from 
the  vicinity  of  airfields.  In  a  series  of  field  experiments  (Boucher, 
1960)  it  was  shown  that  warm  fogs  could  be  acoustically  treated  to 
increase  visibility  by  as  much  as  100  percent.  This  success  seemed  to 
be  limited  to  quiescent  conditions  and  no  improvements  in  visibility 
were  noted  when  wind  speeds  were  greater  than  around  5  to  10  m.p.h. 
Neverthe] ess ,  the  important  point  was  the  demonstration  that  cloud 
droplet  spectra  can  be  modified  by  the  application  of  acoustic  energy. 
This  fact  suggests  the  mechanism  of  acoustic  agglomeration  has  possible 
ramifications  for  the  field  of  meteorology  in  the  effects  of  thunder  on 
cloud  processes  and  as  a  potential  tool  for  use  in  weather  modification 
efforts . 
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The  purpose  of  this  thesis  is  to  briefly  revi?t,  m'.tirt'. 
acoustical  theory  au(?  aerosol  mechanics  relating  to  acoustit 
agglomeration  end  to  describe  the  results  of  t-wo  ijxuj\e'-“icsi  model' 
simulating  acoustically  stimulated  cloud  droplets. 
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Giftpter  II 

THE  SOUND  FIELD  AND  THE  SINGLE  AEROSOL 

The  theory  of  acoustic  agglomeration  has  been  reviewed  by  Mednikov 
(1965),  Shirokovs.  (1970),  and  Shaw  (1978).  In  this  chapter  the  theory 
is  outlined  in  some  detail  since  it  is  a  topic  not  frequently 
encountered  in  meteorology.  The  discussioi;  closely  follows  Mednikov 
(1965).  The  theory  discussed  here  is  limited  to  vhe  traveling  wave 
phenomena  since  in  the  free  atmosphere  significant  standing  sound  wave 
phenomena  are  not  expected  to  occur.  In  the  following  pages  the  theory 
is  developed  by  considering  the  structure  of  the  sound  field  and  then 
examining  the  forces  a  single  aerosol  in  that  field  encounters.  It 
should  be  noted  that  the  word  aerosol  will  be  used  generically  to  refer 
to  any  dispersed  matter  suspended  in  air  and,  in  particular,  is  used 
synonymously  with  droplet  throughout  this  thesis. 

A  discussion  of  the  interactions  between  droplets  in  a  sound  field 
is  deferred  until  the  next  chapter, 

2.1  IM  SOUND  FIELD 

The  coagulation  or  agglomeration  of  aerosol  particles  is  described 
by  Shirokova  (1970)  as  a  series  of  basic  subprocesses,  r,amely,  particle 
convergence,  collision  of  the  particles,  coalescence  of  the  particles, 
and  precipitation  of  the  agglomerates.  Here  we  are  concerned  with  the 
forcing  of  these  subprocesses  due  to  the  application  of  sonic  waves  to 
the  medium. 
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Plane  progressive  soond  waves 


In  an  ideal  gas  when  the  sound  intensity  level  is  less  thi*n  about 
110  db  the  sound  field  aay  be  described  in  terms  of  variations  in  the 
pressure  p',  the  velocity  u',  and  the  fluid  displacement  a'  of  the 
medium.  These  three  periodic  functions  are  expressed  in  the  familiar 
forms 

p '  =  P  coswt  ,  ( 2.1 ) 
a 


and 


n'  =  D  cosiot  3 

a 


(2.2) 


a'  =  A  sinut  (2.3) 

a 

where  P  ,  D  ,  and  A  are  the  pressure,  velocity,  and  displacement 
a  a  a 

amplitudes  respectively,  and  the  frequency  in  rad  s  ^  (u)  or  in  Hz  (f) 
is  related  to  the  wave  speed  c  and  the  wave  number  k  by 

(I)  =  2nf  =  ck  .  (2.4) 


The  speed  of  the  sound  wave  propagation  is  given  by 

c  =  (rRT)^^^  =  (pr/p  )^''^  .  (2.5) 

A 

where  y  is  the  ratio  of  the  specific  heats,  R  is  the  universal  gas 
constant,  T  is  the  mean  temperature  of  the  fluid,  p  is  the  ambient 
pressure  of  the  fluid,  and  is  the  density  of  the  fluid. 

The  y  locity,  pressure,  and  displacement  amplitudes  are  related  to 


the  sound  intensity  by 


23 


i/2 


(2,  )) 


D 

& 


and 


P  =  p  eD 
a  & 


(2.7  I 


2.8} 


where  the  sonnd  intensity  I,  which  is  a  measnre  of  acoustic  energy  i  nx, 
is  given  by  the  expression 


I 


Pa‘= 


(2.  ») 


The  sound  intensity  ir  usually  expressed  as  a  sonnd  intensity  level  i 
decibels  (db)  given  by 


IL  =  10  log^gd/I^) 


(2.10) 


where  =  10  W  m  ^  (10“^^  W  cm~^) . 

This  description  of  a  sonnd  wave  assumes  a  plane  wave  front  in  an 
ideal  medium  and  the  relationships  are  derived  by  linearization  of  the 
equations  governing  fluid  motion  where  second  and  higher  order  terms 
have  been  neglected.  For  situations  in  which  the  intensity  of  the  sound 
is  large  there  are  higher  order  effects  which  must  be  considered.  Some 
of  these  effects  include  inhomogeneities  in  the  sound  field,  distortion 
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in  the  shape  of  the  wave,  increased  absorption  by  the  nedium,  radiatios 
pressure,  acoustic  wind,  and  acoustic  turbulence.  In  subsequent 
sections  some  higher  order  effects  which  relate  to  the  behavior  of 

aerosol  particles  in  the  sound  field  will  be  discussed. 

2.2  TIE  MOTION  ^  ^  ARROSOT.  IN  A  SOUND  FTET.H 

In  this  section  the  forces  acting  on  a  single  droplet  are 

considered  as  an  introduction  to  the  effects  of  sound  on  aerosol 

mixtures.  The  forces  acting  on  the  isolated  aerosol  can  be  classified 

into  two  categories,  namely,  viscous  forces  and  drift  forces.  The 
viscous  effects  of  the  fluid  produce  entrainment  of  the  aerosols  in  the 
fluid  motion  while  the  drift  forces,  which  arise  due  to  second  and 
higher  order  effects  of  scnnd  in  the  fluid,  produce  a  steady  drift  of 
the  aerosols.  Both  of  these  categories  are  discussed  in  the  following 
pages . 


2.2.1  Entrainment 

Entrainment  refers  to  the  degree  to  which  aerosol  particles 
participate  in  the  vibrational  motion  of  the  fluid.  When  a  sound  field 
is  set  up  in  a  volume  of  air  the  fluid  particles  will  begin  to  oscillate 
in  a  periodic  manner.  Aerosols  contained  in  the  volume  will  also  tend 
to  oscillate  but  will  not  be  fully  entrained  in  the  oscillatory  motion. 
The  amount  of  entrainment  may  be  approximated  from  the  Stokes  equation 


for  the  motion  of  a  particle  given  by 


(2.13  ) 


a  —  =  67iT|a(n' 
°  dt 


) 


where  t.  is  the  aerosol  speed,  m,  is  the  aerosol  mass,  tj  is  the 
d  u 

dynamic  viscosity  coefficient  of  the  flnid,  and  s  is  the  radius  of  the 
aerosol.  The  relaxation  time  for  the  droplet,  i.e.,  the  time  required 
for  the  droplet  and  the  <iir  to  come  to  equilibrium  if  either  is 
impulsively  pushed  past  the  other,  is  defined  by 


T 


2p 


d  2 
a 


(2.12) 


Since 


m  =  -p  na 
3 


equation  (2.11)  may  be  rewritten  using  (2.2)  and  (2.12)  to  obtain 


‘‘“•1  „  . 

r —  +  u,  =  U  sinoit 

,  d  a 

dt 


(2.13) 


which  has  a  general  solution 


n  sin(a>t-E)  U  wt 

^  _ _ _  _ K  _  -t/i 

^  ^  r  2  2,i72  ,  ^  "2~2  ® 

(i  +  t!)r)  l+u)r 


(2.14) 


wher’  5  is  a  phase  shift  such  that 


tan^  =  wT 


(2.15) 
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'rhis  p^css  sixift  arises  due  to  tiie  inertia  oi  the  droplet  which  prevents 


it  from  moving  exactly  vita  the  air  motion.  The  second  R.H.S.  term  of 
(2.14)  approaches  zero  after  a  time  t  =  T  so  that  the  droplet  motion 
mav  be  described  by 


fj  sin(ut~5} 


(1  + 


2“2,172 

W  T  ) 


(2.16) 


The  droplet  and  the  air  each  reach  maximum  velocity  when  the  cos  term  in 
(2.2)  and  sin  term  in  (2.16)  each  equal  one.  The  maziranm  speed  of  the 
drop  is  then 


u 


d 


U 

a 


(1 


■  TTili 

+  (J)  T  ) 


Equation  (2.16)  may  be  rewritten  as 


sin((i)t-f) 

d  d  a 

where 


(2.17) 


(1  + 


2  2,1/2 

(1)  X  ) 


(2.18) 


Equation  (2.18)  may  be  rewritten  n.*-ing  equation  (2.15)  to  get 

=  ccs^  .  (2.18a) 

Substituting  for  oj  and  r  from  (2.4)  and  (2,12)  yields  the  two 
equations  for  entraiiunent 


1 


P  j  ®  f  T  -j  /  ^ 

[  1  - 

9t\ 


and 


tan$  = 


4T?p^a^f 

9n 


(2o20) 


Thns  the  degree  of  entreinment  is  determined  by  the  size  of  the  term 

2  2  2 
(4np^8  f/9-n)  and  in  particular  by  the  product  a  f. 

Another  useful  parameter  is  the  flow  around  factor,  i.e,,  the  flow 
of  the  medium  relative  to  the  entrained  droplet.  An  expression  may  be 
derived  by  subtracting  the  drop  speea  (2,17)  from  the  speed  of  the  air 
(2.2) 


u  ,  =  U  sinu't  -  u.U  sin(u)t-£)  . 
ad  a  a  ^ 

Now  using  trigonometric  identities  and  equation  (2.18a)  the  flow  around 
speed  is  written 

2 

u  ,  =  (siniut  -  sinwteos  E  -r  cosiots inScosE) U 
ad  ^  ^  ^  a 

or 

u  ,  =  p  U  cos(u)t~E)  (2.21) 

ad  '^a  a 

where 


tax 


=  sinE  = 


2  2,1/; 


(1  +  (t)  r  ) 


Additionally,  it  can  be  seen  from  (2.18a)  and  (2,22)  that 


(2.22) 
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(1  +  0)  T  )■'' 


Aaiit  ional  ly,  it.  can  be  seen  from  (2.18a)  and  (2.22)  that 


2  .  2  _  , 
u,  =  1  . 

a 


(2.23) 


Fig.  1  shows  the  plots  of  entrainment  and  flow  aronnd  for  various 
droplet  sizes  and  sound  frequencies.  It  can  be  seen  from  the  figure 
that  for  small  droplets  at  low  sound  frequencies  entrainment  is  almost 
complete  while  for  large  droplets  at  high  frequencies  the  amount  of 


entrainment  is  minimal. 
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Figure  1:  Entrainment  and  Flow-aronnd  Factors  for  Various  Sound 

Frequencies.  (Solid  -  entrainment,  dashed  -  flow-around; 
after  Shaw  and  Rajendran,  1979.) 


2.2.2  Particle  drift 

Aerosol  particle  drift  may  play  an  important  ro)e  in  acoustic 
agglomeration  by  moving  aerosols  close  to  one  another  so  that  mechanisms 
discussed  in  Chapter  III  can  cause  collisions  to  occur.  Among  the 
factors  contributing  to  particle  drift  in  a  trc’/eling  wave  are  radiation 
pressure,  periodic  viscosity  changes,  vibrational  phase  differences,  end 
distortion  of  the  sound  wave.  These  will  be  described  in  turn. 

When  8  sound  wave  encounters  any  .bstacle  in  the  fluid  'e  wave 
energy  is  scattered.  In  the  case  of  a  drop  some  of  the  energy  is 


scattered  in  the  direction  of  propagation.  There  is  a  momentum  transfer 
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associated  #ltJi  this  process  and  since  more  momentum  is  transferred  to 
th'*  drop  by  the  trave  than  is  transferred  by  the  drop  in  the  direction  of 
propagation,  there  remains  a  net  force  on  the  drop  which  is  the  momentum 
bacL-scattered  toward  the  source.  Expressions  for  the  radiation 
pressure  force  on  a  drop  have  been  derived  by  King  vl934)  and  Landau  and 
Lifshitz  (1959).  The  radiation  pressure  force  for  spherical  particles 
is  given  by 


R 


11  0) 
— Jt(-) 

9  c 


4  6  2„ 
a  u  E 
'^a 


(2.24) 


where  E  =  2I/c  is  the  sound  energy  density. 

The  next  two  drift  forces  were  proposed  by  Westervelt  (1950)  in  a 
paper  dealing  with  the  mean  pressure  and  particle  velocity  in  a  sound 
wave.  Tn  that  paper  Westervelt  treated  the  two  forces  as  a  single  force 
wnich  he  termed  a  Stokes  force.  Mednikov  (1965)  separated  the  two 
forces  attributing  one  to  changes  in  viscosity  and  the  other  to  phase 
differences  in  the  vibration  of  tic  fluid  particles. 

First  the  drift  force  due  to  changes  in  viscosity  is  discussed. 
Because  of  the  periodic  adiabatic  compressions  and  rarefactions  of  a 
sound  wave  there  occur  periodic  temperature  changes.  Therefore,  since 
it  is  a  function  of  temperature,  viscosity  increases  during  compression 
and  decreases  during  rarefaction.  During  compression  a  droplet  is  being 
moved  in  the  direction  of  propagation  and  because  the  viscosity  is 
greater  the  droplet  is  more  fully  entrained  in  the  fluid  motion.  The 
exact  opposite  process  occurs  during  a  rarefaction.  ihe  result  is  a  net 
force  in  the  direction  of  propagation  given  by 
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(2,25) 
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.?n(Y-l)-a^2g 


Thfe  drift  force  which  arises  due  to  phase  differences  in  the 
vibration  of  adjacent  fluid  particles  is  very  complex  end  difficult  to 
explain  in  either  a  dynamic  or  kinematic  sense.  For  a  detailed 
discussion  the  reader  is  referred  to  Westervelt  (1950),  Westervel^ 
(1951),  and  Mednikov  (1965).  An  expression  for  the  resulting  force  on 
aerosols  suspended  in  the  fluid  is  given  by 


2 

F  =  -6n-aii  E  . 
P  c  » 


(2.26) 


The  last  drift  force  discussed  is  the  one  due  to  distortion  or  wave 
asymmetry.  High  intensity  sound  waves  can  become  disterte.-  and  assume  a 
sawtooth  structure.  During  a  compression  the  fluid  particles  are  moving 
in  the  direction  of  propagation  and  thus  pick  up  some  of  the  wave  speed. 
During  rarefaction  the  particles  move  opposite  to  propagation  with  the 
pressure  induced  velocity  less  some  of  the  wave  velocity.  As  a  result 
the  wave  slope  is  usually  steeper  during  the  compression  portion  of  the 
wave  than  during  the  rarefaction  portion,  and,  therefore,  the  velocity 
usually  increases  rapidly  in  the  direction  of  propagation  while 
increasing  more  slowly  in  the  opposite  direction.  Since  an  aerosol  in 
the  regior.  of  rapidly  increasing  speed  will  be  less  fully  entrained  than 
when  it  is  in  the  region  of  slowly  increasing  speed  there  results  a  net 
force  in  the  direction  opposite  to  propagation.  This  force  may  be 
expressed  as 
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(2.27) 


F.  =  -61i.,siiiX  a  u  E  . 
h  2  'a 

rhere  relative  velocity  amplitude  of  the  second  harmonic  and 

X  is  the  phase  angle  of  the  second  harmonic. 

The  relative  importance  of  the  drift  forces  under  various 

conditions  can  be  seen  by  evaluating  (2.21)  through  (2,25)  using  typical 

2 

values  for  the  variables  and  factoring  the  common  term  afi  E  from  each 

a 

expression.  Table  1  shows  relative  values  (cm)  of  the  scaled  forces  for 

4-1  -12 

three  droplet  sizes.  Here  f  =  1000  Hz,  c  =  3x10  cm  s  ,  v  -  10  cm 
s  and  Y  =  1.4  .  The  values  for  "^2  sinX  are  taken  to  be 

unity  in  order  to  obtain  the  maximum  value  of  the  expression.  It  is 
clear  that  the  largest  force  is  that  due  to  as3Tnmetry  of  the  sound  wave 
when  the  sound  wave  is  significantly  distorted.  When  distortion  is  not 
a  factor  then  the  drift  due  to  particle  vibration  phase  differences  is 
generally  the  largest.  Also  it  is  clear  that  for  applications  of 
interest  in  this  study  the  radiation  pressure  force  is  negligble. 
However,  being  strongly  dependent  on  particle  radius,  this  force  can 
become  dominant  for  large  particles  and  is  responsible  for  the 
phenomenon  of  suspended  objects  in  a  sound  field  (St.  Clair,  1949). 
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TABLE  1 


Relative  Values  of  Drift  Forces 
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9  c 

-23 

7.0x10 

7.0x10“^* 

-13 

7.0x10 

Viscosity 
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Phase 

V 

V 
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c 
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Asymmetry 

-dh^a  sinX 
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-5.4xl0~^ 
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2.3  A  SINGLE  DROP  MODEL 

A  single  droplet  numerical  model  is  used  to  examine  the  motions  of 
a  drop  within  a  sound  field.  As  noted  earlier  this  model  simplifies  the 
problem  of  Jropiet  behavior  in  an  acoustic  field  since  it  neglects  the 
presence  of  other  droplets  or  particles  which  certainly  force  motions 
upon  the  droplet  under  consideration.  But,  this  deficiency  aside,  it  is 
of  interest  to  study  the  motions  of  differen.  sized  drops  under  the 
influence  of  sound. 

2.3.1  The  Model  Equations 

The  sound  intensity  required  for  the  agglomeration  of  aerosols  is 
high  enough  that  distortion  of  the  sound  wave  is  significant.  Thus,  for 
the  purposes  oi'  this  example,  the  radiation  pressure,  viscosity  change. 
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and  vibrational  phase  forces  are  neglected  based  on  the  results 
described  in  Table  1. 

In  the  previous  section  it  was  shown  that  the  motion  of  a  droplet 
due  to  the  vibration  of  the  air  is  given  by  equation  (2.16).  To 
describe  the  motion  due  to  the  drift  forces  it  is  assumed  the  drop 
relaxation  times  are  very  small  while  changes  in  the  drift  forces  are 
much  longer.  These  assumptions  approximate  an  equilibrium  where  the 
drift  forces  are  in  a  continuous  balance  with  the  Stokes  resistive 
force.  Thus  the  velocity  due  to  each  force  is  obtained  by  equating  the 
force  to  the  Stokes  drag  force 


F 

ST 


=  6nT]au^ 


Then  the  approximate  total  velocity  of  the  drop  is  given  by 


u,  =  U  sin(kx-wt) 
d  a 


h^sinX  2 

- ail  E 

nn 


(2.28) 


where  it  is  assumed  that  the  phase  shift  of  the  second  harmonic  is  n/2. 

In  order  to  examine  the  validity  of  the  Stokes  flow  assumption  a 
Reynolds  number  for  a  drop  in  an  acoustic  field  is  defined  b;-  the 
expression 

2aU 

Re^^  =  - -  (2.29^ 

V 


where  v  is  the  kinematic  viscosity  coefficient  and  is  the  maximum 
velocity  of  the  droplet  given  by  The  Reynolds  numbers  for 
drops  at  different  sound  intensity  levels  are  plotted  in  Fig,  2.  The 
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solid  curve  ic  the  figure  is  a  plot  for  drops  falling  at  their  terEinal 
velocities.  The  dashed  curves,  by  increasing  slope,  represent  sound 
intensity  levels  of  100  dh,  120  db,  140  db,  and  160  db,  respectively. 
As  an  example,  the  Reynolds  number  is  less  than  one  for  drops  up  to  10  p 
radius  at  a  sound  intensity  level  of  140  dh.  For  situations  in  which 
the  Reynolds  number  is  greater  than  one,  an  Osees  correction  (Mednikov, 
1965)  is  applied  to  the  Stokes  velocity  hy  use  of  the  equation 
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(2.30) 
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2.5 Pesalts 

A  Eodel  based  oa  eqaation  (2,28)  has  been  used  to  esaisiae  the 
iBOtioas  of  various  sized  droplets.  In  the  eransples  discussed  here  the 
values  of  the  atmospheric  and  sound  variables  included  p  =  765  mb,  T  = 
283  K,  IL  =  140db.  and  f  =  1000  Hz, 

In  Fig.  3  through  Fig.  5  plots  of  displacement  vs.  time  are  shown 
for  drops  of  1,  10,  and  20  p  respectively.  The  scale  of  tii*>  plots  is 
the  same  in  each  figure.  It  is  evident  that  the  greatest  amplitude  of 
vibration  occurs  in  the  case  of  the  1  p  droplet  while  there  is  no 
apparent  drift  for  this  same  droplet.  Further,  it  is  clear  from  these 
figures  that  as  drop  size  increases,  vibrational  amplitude  decreases  and 
drift  velocity  increases. 

In  this  section  the  motions  of  isolated  droplets  were  examined.  In 
the  next  chapter  the  discussion  focuses  on  the  more  important  case  where 
adjacent  aerosols  induce  motions  on  one  another. 
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Fignre  3:  Motion  of  a  1  micron  droplet 


Chaptex  III 


INI^RACTIONS  OF  AEROSOLS  IN  A  SOUND  FIELD 

The  motions  discnssed  in  the  previotjs  chapter  do  not  take  into 
account  the  presence  of  nearby  aerosols  which  can  have  important  and 
even  dominating  effects  on  aerosol  motion.  In  this  chapter  some  of  the 
motions  due  to  drop  interactions  in  an  acoustic  field  are  discussed. 

Mednikov  (1965)  classifies  four  types  of  motions  due  to  interaction 
as  orthokinetic,  parakinetic,  attractional ,  and  pulsational.  In  the 
following  description  the  interactions  are  more  generally  classified  as 
orthokinetic,  hydrodynamic,  and  turbulent  where  the  parakinetic  and 
attractional  interactions  are  grouped  under  the  hydrodynamic 
classification.  Tbe  pulsational  interaction  will  be  called  turbulent 
interaction  in  this  thesis. 

3.1  ORTHOKINETIC  IbriERACTION 

The  orthokinetic  interaction  is  that  which  occurs  between  drops 
whose  centers  lie  along  a  line  which  is  parallel  to  the  direction  of 
motion.  This  motion  may  be  the  result  of  drift  or  some  turbulent 
pulsation  but  the  most  significant  contributor  is  the  vibrational  motion 
induced  by  the  periodic  sound  wave.  An  expression  for  the  relative 
velocity  of  two  droplets  may  be  found  by  subtracting  the  respective 
velocities  given  by  (2.17)  to  get 
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To  find  the  maximnm  difference  in  speed  the  above  equation  is 
differentiated  with  respect  to  time  and  set  equal  to  zero  to  get  an 
expression  for  the  frequency.  Using  equations  (2.18a)  and  (2.22)  the 
result  is  written 


tanutp 


2 

fi2. 


■  »^d2^2 


where  t^  is  the  time  at  which  the  maximum  difference  occurs.  Now  the 
maximum  relative  velocity  i?-  termed  U*  and  is  given  by 

or 


U*  -  [^^^(sinwtQCOs^^  -  cosojtQsin^j^) 

~  ~  coswtgsin^^)  ^ 


Using  (2.18a)  and  (2.22)  gives 


U-*  =  [sina,tQ(^2^-pJ^)  -  cosa)tQ(n^^p^2■^^^2^2^^  ^a 


Thus , 


U*  = 


'*12''. 


where 


2  2. -1/2., 
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2  2, -1/2 

“  ^2^ 


(3.3) 


Equations  (2.18)  and  (2.22)  have  been  nsed  to  rewrite  |ij^2  terici  of 
(D  and  %.  By  taking  the  derivative  of  (3.1)  with  respect  to  w  and 
setting  it  equal  to  zero,  an  optiminn  frequency  for  agglomerating  two 
droplets  by  the  orthokinetic  mechanism  is  found  to  be 


^OPT 


9il 


(3.2) 


While  giving  useful  information  about  the  frequency  to  use  on  an  initial 
spectrum,  this  expression  is  of  limited  practical  applicability  because 
the  droplet  sizes  begin  changing  as  soon  as  coagulation  begins. 


3.2  HYDRODYNAMIC  INTERACTIONS 

The  physical  basis  for  the  hydrodynamic  (attractional  and 
parakinetic)  interactions  is  well  known  and  has  been  described 
theoretically  by  several  authors  (  e.g.  Lamb,  1932).  It  consists  of  the 
motions  introduced  as  a  result  of  the  superposition  of  the  flow  fields 
around  each  droplet.  In  the  simplest  case  there  are  two  spheres  of  the 
same  dimensions  located  in  a  fluid  and  separated  by  some  distance.  If 
the  fluid  flows  past  the  spheres  such  that  the  streamlines  are  parallel 
to  the  line  joining  the  centers  of  the  spheres,  the  spheres  will  move 
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toward  one  another.  Thus  the  attractional  interaction  in  a  sound  field 
pertains  to  the  approach  and  collision  of  droplets  whose  centers  are 
essentially  aligned  with  the  streamlines  of  the  acoustically  induced 
fluid  flow. 

The  parakinetic  interaction,  on  the  other  hand,  pertains  to  motion 
that  occurs  between  two  droplets  whose  centers  are  along  a  line  not 
parallel  to  the  streamlines.  The  idealized  motion  as  described  by 
Mednikov  (1965)  is  a  zig-zagging  of  the  smaller  droplets  in  the  vicinity 
of  a  larger  drop.  The  zig-zagging  motion  is  such  that  droplets  to  the 
front  or  back  of  the  collector  tend  to  become  centered  along  the 
collector's  line  of  motion.  Those  droplets  which  are  initially  along  a 
flank  of  the  collector  tend  to  be  moved  transversly  away  from  the 
collector  until  they  are  caught  up  in  streaming  of  the  fluid  past  the 
collector  and  are  carried  forward  or  aft  of  the  collector  where  they  may 
begin  the  centering  motions  described  previously. 

Since  the  hydrodynamic  interactions  are  important  to  the  modeling 
of  the  agglomeration  process,  an  expression  for  the  speed  of  approach  of 
two  droplets  in  an  zcoustic  field  is  needed.  The  agglomeration  model 
that  will  be  discussed  in  Chapter  TV  is  based  on  the  work  of  Chou,  et. 
al.,  (1981)  who  used  an  expression  derived  from  theory  by  Dianov,  et. 
al.,  (1968).  This  is  a  general  expression  for  the  relative  velocity  of 


two  spheres,  not  necessarily  of  the  same  size,  under  the  conditions  of 
Oseen  flow  in  an  acoustic  field.  The  relative  velocity  is  given  by 
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where  s  is  the  separation  distance  of  the  droplets.  Other  terms  in 


(3.3)  are  defined  by 
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5  =  tan"^( - 5 - -2  ) 


where  ji  and  (i,  are  given  by  eijnations  (2.18)  and  (2.22)  respectively. 
The  first  three  terms  of  (3.3)  represent  the  ve.lccity  due  to  the  mutual 
distortion  of  the  flow  fields  around  the  two  droplets.  The  fourth  term 
represents  the  velocity  due  to  the  vibrational  phase  difference  of  the 
droplets . 

Figure  6  shows  plots  of  the  relative  velocities  of  various  droplet 
pairs  due  to  gravitational  sedimentation  and  the  acoustic  hydrodynamic 
interaction.  The  relative  velocities  due  to  gravity  were  computed  by 
subtracting  the  terminal  velocities  of  the  vrc*'lets. 


relative  velocities  were  computed  from  equation  (3.3).  For  this 
acoustic  relative  velocity  calculation  the  droplet  separation  is  100  u, 
the  sound  intensity  level  is  140  db ,  and  the  frequency  is  1000  Hz.  The 
top  pair  of  solid  and  dashed  curves  is  for  a  droplet  with  a  radius  of  22 


1  f  otiH  1  0 


respect  ivelv.  T1i:j  relative 


velocity  of  two  droplets  falling  at  their  terminal  velocities  approaches 
zero  the  closer  the  droplets  are  to  the  same  size.  However,  in  an 
acoustic  field  the  relative  velocity  increases  as  the  ratio  of  the  drop 
sizes  approaches  unity.  As  can  be  seen  in  the  figure,  larger  droplets 
may  be  expected  to  approach  one  another  in  an  acoustic  field  at  a  rate 
significantly  greater  than  experienced  under  the  force  of  gravity  alone. 
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HYDRODYKAMIC  AND  TERICNAL  REUTtVE  VELOCmES 
AS  A  FUNCTION  OF  RAORTS 


HfiDius  tM:c=CNs: 


Figure  6:  P^lative  Velocities  Dae  to  Gra-vity  and  Acoustics.  (Solid 
’urve  is  for  drops  in  an  acoustic  field,  d^sLed  curves  fc 
drops  falling  at  terminal  velocity.) 


% 


3,3  TURBULEOT  INTERACTIONS 


The  onset  of  significant  acoustic  turbulence  does  not  occur  until 
sound  intensities  reach  a  relatively  high  level.  With  respect  to  the 
agglomeration  process  the  onset  of  the  turbulence  effect  does  not  occur 
until  the  sound  intensity  reaches  levels  around  160  db  (Chou,  et.al., 
1981).  Since  sound  intensities  that  high  are  generally  not  expected  to 
be  present  over  large  cloud  volumes,  the  models  discussed  in  this  thesis 
do  not  incorporate  any  turbulent  processes.  However,  a  short  discussion 
of  the  turbulent  effects  is  provided  for  general  background  in  the 
subject  of  acoustic  agglomeration. 

Thv.-  turbulent  or  pulsational  interactions  are  due  to  motions 
arising  from  the  acoustically  induced  turbulence  in  the  flow  field. 
There  are  essentially  two  interactions.  The  first  is  described  as 
inertial  and  is  characterized  by  a  turbulence  length  scale  small 
compared  to  the  aerosols.  Thus,  the  droplets  are  not  fully  entrained  in 
the  fluid  motion  and,  since  the  motion  field  itself  is  inhomogeneous, 
the  droplets  simply  cannot  follow  the  rapid  fluctuations  in  the  fluid 
motion  and  collisions  occur  in  a  stochastic  manner. 

The  second  turbulent  interaction  is  a  diffusion  process  in  which 
the  length  scale  of  the  turbulence  is  large  compared  to  the  aerosols. 
Her.;  droplets  become  fully  entrained  in  the  fluid  motion,  but,  adjacent 
droplets  have  collisions  because  they  have  differing  relative  velocities 
due  to  the  spatial  inhomogeneities  of  the  turbulent  motion. 
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3.4  SIBffiAKY 


Thongh  all  the  motions  and  interactions  described  to  this  point 
certainly  play  some  part  in  the  process  of  aeonstic  agglomeration  a 
complete  theory  of  the  process  has  not  been  fnlly  developed.  It  was 
first  believed  (Brandt  and  Hiederoann,  1936)  t^tt  the  orthokinetic  effect 
▼as  the  cause  of  agglomeration  but  since  that  effect  requires  size 
differences  among  the  particles  and  agglomeration  occurs  in  monodisperse 
populations  as  well  as  polydisperse ,  the  theory  required  extension. 
Mednikov  (1965)  proposed  that  the  most  significant  cause  of 
agglomeration  was  the  orthokinetic  interaction  while  al<  the  other 
motions  acted  to  refill  the  collision  volume.  Shirokova  (1S7U)  argued 
that  the  orthokinetic  mechanism  is  not  as  important  ss  the  hydrodynamic 
interaction  of  the  particles.  However,  Cheng,  et.  al . ,  (1982) 
demonstrated  that  the  orthokinetic  effect  was  dominant  for  the 
agglomeration  of  submicron  particles  when  there  was  no  acoustic 
turbulence  present.  In  his  review  article  Shaw  (1978)  states  that  the 
hydrodynamic  interaction  may  be  the  most  important  cause  of 
agglomeration  particularly  for  relatively  large  sized  particles.  Shaw 
and  Tu  (1979)  demonstrated  experimentally  that  the  hydrodynamic 
mechanism  is  important  for  ap..osols  greater  than  about  0.5  p  in  radius. 

In  t>'e  next  chapter  a  model  of  the  agglomeratio ,  process  is 
described  and  the  results  of  numerical  experiments  for  atmospheric 
conditions  are  discussed. 
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Chapter  IV 


MODELING  CHANGES  IN  CLOUD  DROPLET  CONCENTRATIONS 

Recent  work  on  acoustic  agglomeration  in  the  United  States  has  been 
done  by  Shaw  and  his  associates  at  the  Laboratory  for  Power  and 
Environmental  Studies  at  the  State  University  of  New  York  at  Buffalo. 
They  have  developed  numerical  models  which  predict  changes  in  aerosol 
spectra  in  good  a^^eement  with  their  experimental  observations.  The 
work  to  be  described  in  this  chapter  is  based  on  the  work  of  Chou,  Lee, 
and  Shaw  (1981)  in  which  the  results  of  a  numerical  model  designed  to 
simulate  the  acoustic  agglomeration  process  in  an  aerosol  mixture  are 
compared  to  the  results  of  laboratory  experiments.  Their  numerical 
model  was  based  on  the  theory  of  acoustic  agglomeration  and  its  results 
were  in  good  agreement  with  the  experimental  results.  In  the  research 
described  here  a  similar  numerical  model  has  been  used  to  examine  the 
effects  of  acoustic  agglomeration  in  an  environment  more  characteristic 
of  clouds. 


A  •;  'fHF  unTTi'TTnMC 

4.1.1  The  Stochastic  Collection  Equat ion  (SCE) 

The  moael  used  here  is  based  on  the  stochastic  collection  equation 
given  by 
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where  K  is  a  collection  kernel,  i  and  x'  are  droplet  volnmes,  and  n  is 
the  nnnher  concentration  per  drop  volnsae  at  time  t.  The  first  R.H.S. 
term  represents  the  gain  in  nnmber  of  droplets  of  volume  x  dne  to 
collisions  between  droplets  with  volnmes  less  than  that  of  x.  The 
second  R.H.S.  term  represents  a  loss  of  droplets  of  volume  x  dne  to 
collisions  with  drops  of  all  other  volnmes  resulting  in  drops  of  volume 
greater  than  i.  The  physics  of  the  collection  kernel  will  be  discussed 
in  the  next  section. 

4.1.2  The  Collection  Kernel 

Previous  work  on  acoustic  agglomeration  has  been  mainly  concerned 
with  highly  concentrated  aerosols  contained  within  small  volumes  over 
which  it  is  feasible  to  maintain  high  sound  intensities.  In  the  Chou 
study,  which  was  concerned  with  conditions  such  as  those  just-  listed, 
four  collection  kernels  were  used.  There  was  one  kernel  defined  for 
each  of  the  following  acoustic  aerosol  interactions:  orthokinetic , 
hydrodynamic,  diffusional  turbulent,  and  inertial  turbulent.  However, 
practical  limitations  suggest  that  it  is  not  feasible  to  maintain  high 
intensity  level  0150  db)  sound  over  the  volumes  encountered  in  working 
with  atmospheric  clouds.  For  this  reason  the  model  used  in  this  work 
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neglects  the  effects  of  turbulent  processes  which  become  important  at 
sound  intensity  levels  of  around  160  db.  Thus,  the  acoustic 
agglomeration  kernels  used  in  this  work  are  for  the  hydrodyT- -jic  and 
orthokinetic  processes.  In  addition,  a  kernel  for  gravitational 
sedimentation  is  added  to  simulate  the  coalescence  process  that 
continues  after  the  simulated  sound  is  terminated. 

The  concept  of  a  collection  kernel  is  based  on  the  volume  swept  out 
per  unit  time  by  a  colliding  droplet  pair.  The  volume  is  called  a 
collision  volume  and  its  cross  section  is  the  geometric  cross  section 
multiplied  by  a  collection  efficiency.  The  length  or  height  of  the 
volume  is  the  distance  traveled  by  the  approaching  drops  per  unit  time. 
The  most  generally  recognized  collection  kernel  is  that  for 
gravitational  sedimentation  or  geometric  sweep  out  given  by 
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where  is  the  collection  efficiency  due  to  gravity  and  v^  is  the 
terminal  velocity  of  the  droplet. 

The  two  acoustic  collection  kernels  used  here  are  the  same  as  those 
used  in  Chou,  et.  al .  (1981).  The  orthokinetic  collection  kernel  was 
derived  by  Mednikov  (1965).  In  this  case  the  cross  sectional  area  of 
the  collision  volume  is  the  same  as  described  above.  But  the  length  of 
the  volume  is  twice  the  relative  displacment  amplitude  of  the  drops 
multiplied  by  the  frequency  or  number  of  oscillations  per  unit  time. 
The  collection  kernel  may  be  written 
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TfLexe  is  the  relative  amplitude  of  vibration  of  the  two  droplets. 

After  one  complete  vibration  the  large  drop  is  assumed  to  have  kollected 
all  the  little  drops  within  the  agglomeration  volnme.  However, 
additional  small  drops  may  move  into  the  volume  doe  to  any  of  the  drift 
or  interaction  motions  discussed  is  earlier  chapters.  Therefore,  aj 
additional  coefficient,  called  a  refill  factor,  is  necessary  to  account 
for  small  droplets  moving  back  into  the  voinme  after  each  oscillation. 
Thus, 

2  =  2p7iEQg(aj^+a2)  A^^f 

where  p  is  the  refill  factor.  Now,  using  Aj^2  ~  ^12^a  equations 

(2.4)  and  (2.8)  the  orthokinetic  kernel  is  written 
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(4.3) 


where  is  the  number  concentration  of  large  droplets,  N^  is  the  total 
concentration,  and  is  the  relative  entrainment  coefficient  given  by 
equation  (3.1).  In  this  final  form  Giou  has  nsed  p  =  to  define 
the  refill  factor.  The  validity  of  this  choice  is  assumed  to  be 
empirically  based  since  no  analytical  expression  for  the  refill  factor 
has  yet  been  found.  However,  as  mentioned  before,  the  results  of  Chou's 
model  ag’-ee  well  with  experimental  results. 

The  second  acoustic  kernel  is  the  hydrodynamic  kernel  which  is 
given  by 
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where  is  tiie  relative  velocity  given  by  (3.3).  The  physical  basis 
of  this  kernel  is  directly  analogous  to  that  of  the  gravitational 
kernel.  Here  the  relative  motion  of  the  two  droplets  is  due  to  acoustic 
forces  rather  than  gravity. 

It  is  appropriate  to  make  a  few  comments  about  the  collision 
efficiencies  in  the  collection  kernel  equations.  As  noted  above  each  of 
the  collision  efficiencies  is  different  for  each  process.  For  the 
gravitational  kernel  the  collision  efficiency  is  that  used  in  the 
Leighton  and  Rogers  (1974)  version  of  the  Berry  and  Reinhardt  (1974a) 
SCE  model.  That  is  the  Shafri”  and  Neiburger  (1963)  scheme  for  pairs 
with  >  30  p  and  Hocking  and  Jonas  (1970)  for  pairs  with 

In  order  to  find  values  for  the  collision  efficiencies  of  aerosols 
in  an  acoustic  field  Shaw  (Shaw,  1978p  Shaw  and  Rajendran,  1979p  Shaw 
and  Tu,  1979p  and  Cnou,  et.  al.,  1981)  examined  the  dynamical  similarity 
of  aerosols  in  close  proximity  in  a  gravitational  field  versus  aerosols 
in  an  acoustic  field.  He  reasoned  that,  because  of  differences  in 
acoustic  entrainment,  aerosols  in  an  acoustic  field  approach  one  another 
with  some  relative  velocity  in  much  the  same  manner  as  aerosols  falling 
at  their  terminal  velocities.  By  non-dimensional  analysis  of  the 
equations  of  motion  of  aerosols  in  the  two  different  fields,  he  showed 
that  the  collision  process  is  similar  as  long  as  the  actual  radius  of 
the  larger  aerosol  is  replaced  by  an  acoustically  equivalent  radius 
given  by 
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This  formulation  accounts  for  the  fact  that  aerosols  of  a  gives  size 
have  greater  relative  velocities  in  an  intense  acoustic  field  than  due 
to  gravity  alone.  ITius,  by  using  ihe  value  of  a  for  the  large  drop  in 
an  expression  for  the  collision  efficiency  due  to  gravitional 
sedimentation,  the  collision  efficiency  due  to  acoustically  driven 
motion  can  be  found.  In  particular,  Chou  (1981)  uses  the  formula  of 
Langmuir  (1948)  for  the  orthokinetic  collision  eiiiciency  and  the 
formulation  of  Scott  and  Chen  (1970)  for  the  hydrodynamic  collision 
efficiencies.  The  same  formulas  are  used  in  the  model  described  in  this 
thesis . 

It  is  of  some  value  to  examine  the  sensitivity  of  the  acoustic 
collection  kernels  to  variations  in  the  environment.  It  should  be  kept 
in  mind  that  the  discussion  does  not  describe  the  extremely  complex 
agglomeration  process.  The  intent  of  the  following  is  to  construct  some 
framework  for  understanding  the  complete  process  by  examining  very 
limited  interactions.  The  relative  effects  of  environmental  variables 
on  each  kernel  are  illustrated  in  Fig.  7  through  Fig.  13.  These  curves 
were  calculated  for  the  following  environmental  and  sound  field  values:  p 
=  765  mb,  T  =  283  K.  a^  =  12  p,  var  =  1.0,  IL  =  140  db,  f  =  1000  Hz,  and 

-3 

L  =  1.0  g  m  .  Each  of  these  parameters  was  varied  in  turn  to  obtain 
the  kernel  values  depicted  in  the  figures.  The  radii  of  the  droplet 
pair  for  each  kernel  were  selected  arbitrarily  and  the  results  depicted 
vould  change  only  in  magnitude  for  different  droplet  pairs.  In  the 
gravitational  and  orthokinetic  kernels  a^  =  4  p  and  a2  =  20  p  for  all 
cases.  In  the  calculation  of  the  hydrodynamic  kernels  a^  =  a^  and  a2  = 
0.7a^.  In  all  of  the  diagrams  the  solid  line  represents  the 
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ortliokinetic  kernel  and  the  dashed  line  represents  the  hydrodynamic 

kernel.  For  reference  purposes  the  value  of  the  gravitational  kernel 

~6  3  ~1 

for  a^  =  4  p  and  a2  =  20  p  is  =  6.14x10  cm  s  and  is  relatively 
insensitive  to  the  environmental  variations  studied  here. 

Fig.  7  shows  that  both  the  orthokinetic  and  hydrodynamic  kernels 
increase  with  increasing  sound  intensity.  Both  of  these  effects  can  be 
attributed  to  the  increase  in  the  amplitude  and  velocity  of  the  medium's 
oscillations  as  shown  in  equations  (2.8)  and  (2.10).  In  the 
orthokinetic  interaction  the  displacement  of  each  aerosol  increases 
thereby  increasing  the  size  of  the  collision  volume.  Also,  the  relative 
velocities  of  the  collision  pairs  are  greater.  For  the  hydrodynamic 
interaction  the  relative  velocity  of  the  pair  increases  due  to  the 
increased  speed  of  the  fluid  past  the  aerosols. 

The  sensitivity  to  sound  frequency  changes  is  shown  in  Fig  8.  For 
the  orthokinetic  effect  the  kernel  increases  to  a  peak  around  200  Hz  and 
falls  off  asymptotically  from  there.  The  response  of  the  orthokinetic 
kernel  to  changes  in  frequency  is  explained  in  terms  of  the  entrainment 
of  the  drops  which  is  governed  by  equation  (2.19).  For  example,  the 
decrease  in  magnitude  of  the  kernel  with  increasing  frequency  arises 
because  both  droplets  are  less  fully  entrained  in  the  air  motion.  The 
increasing  followed  by  decreasing  b;havior  of  the  curve  occurs  since 
there  are  also  relative  entrainment  effects.  The  maximum  value  for  a 
droplet  pair  occurs  at  the  frequency  given  by  equation  (2.27).  Dianov, 
et.  al.,  (1968)  show  how  the  hydrodynamic  relative  velocity  increases 
asymptotically  with  increasing  frequency.  Since  the  hydrodynamic  kernel 
is  proportional  to  the  relative  velocity  it  increases  in  the  same 
manner , 
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Is  Fig.  9  it  can  be  seen  that  Loth  aconstic  kernels  decrease 
increasing  pressure  and  in  Fig,  10  both  kernels  are  seen  to  increase 
slightly  with  increasing  temperature.  The  effects  of  temperature  and 
pressure  changes  on  both  the  hydrodynamic  and  orthokinetic  processes  are 
a  complex  combination  of  viscosity  and  sound  wave  properties.  For  the 
range  of  atmospheric  values  these  effects  are  not  very  significant. 

In  addition  to  the  sound  and  environmental  variables  discussed  so 
far,  the  acoustic  collection  kernels  are  also  sensitive  variations  in 
the  aerosol  spectrum.  In  the  orthokinetic  case  the  kernel  is  a  function 
of  a  concentration  term  as  can  be  seen  in  equation  (4.3).  Recall  this 
term  parameterizes  the  refilling  of  the  agglomeration  volume  with  small 
droplets.  The  hydrodynamic  kernel  given  by  equation  (4.4)  is  a  function 
of  the  relative  velacity  which  is  a  function  of  droplet  separation  as 
shown  in  equation  (3.3).  Average  droplet  seperation  is  a  function  of 
droplet  concentration.  Thus,  both  acoustic  kernels  are  dependent  on  the 
aerosol  concentration. 

The  aerosol  concentration  is  dependent  upon  the  specification  of 

the  intial  aerosol  spectrum.  In  this  research  work  the  spectra  used  are 

all  of  a  gamma  or  Pearson  Type  III  distribution  and  are  specified  by 

liquid  water  content  L,  initial  mean  radius  a^,  and  shape  factor  or 
2 

variance  o  ,  In  Fig.  11  through  Fig.  13  the  effects  of  changes  in  the 
parameters  that  specify  the-  initial  aerosol  spectrum  are  illustrated. 
In  Fig,  1.1  it  is  shown  that  when  the  initial  mean  radius  is  increased 
the  orthokinetic  effect  increases  and  reaches  a  maximum  around  15  p  and 
then  falls  off  toward  larger  sizes.  In  terms  of  equation  (4.3)  this 
behavior  results  from  the  relative  concentration  term  reaching  a  maximum 


at  =  15  n.  Here  the  hydrodynamic  kernel  behavior  is  an  artifact  of 
the  method  of  calculation,  i.e.,  the  radius  a^  in  equation  (4.4)  is 
equal  to  a^.  In  general,  for  a  given  droplet  pair,  the  hydrodynamic 
kernel  will  decrease  with  increasing  initial  mean  radius  because  the 
droplet  concentration  will  decrease  thereby  increasing  the  mean 
separation  of  the  droplets. 

When  the  initial  spectrum  variance  is  increased  the  orthokinetic 
effect  increases  as  seen  in  Fig  12.  Physically,  increasing  the  variance 
increases  the  relative  numbers  of  large  collectors  and  small  collectible 
droplets.  For  the  orthokinetic  interaction  there  are  more  droplets  with 
different  relative  entrainment  in  the  fluid  motion  resulting  in  more 
collisions.  Again,  in  terms  of  equation  (4.3),  the  relative 
concentration  term  is  increased.  The  hydrodynamic  effect  shows  no 
significant  change.  In  terms  of  equation  (4.4)  any  changes  are  due  to 
variations  in  the  relative  velocity  which  is  inversely  dependent  on  the 
separation  of  the  droplets. 

The  effect  of  increasing  the  liquid  water  content  is  demonstrated 
in  Fig  13.  This  effect  is  a  droplet  concentration  effect  since 
increasing  the  liquid  water  proportionately  increases  the  number  of 
droplets.  The  orthokinetic  kernel  varies  directly  with  the  r  lative 
concentration.  Once  again  the  concentration  effect  on  °.verage 
separation  changes  the  relative  velocity  and  thus  the  hydrodynamic 
kernel . 
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ACOUSTJC  KERNELS  AS  A  EUNCTION  OF 
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Figure  8:  Kernel'  as  a  Function  of  Sound  Frequency.  (Solid  curve  is  for 
orthokinetic ,  aashed  for  hydrodynamic.) 


ACOUSTIC  KERNELS  AS  A  EUMCnON  Of 
AiffiUrriT  PRESSURE 
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Figure  9:  Kernels  as  a  Function  of  Ambient  Pressure.  (Solid  curve  is 
for  ortbokinetic ,  dashed  for  hydrodynamic.) 


ACOUSTIC  KERNELS  AS  A  >’UNCnON  OF 
AMBIENT  TEMPERATURE 


Figure  10:  Kernels  as  a  Function  of  Ambient  lemperature.  (Solid  curv ' 
is  for  orthokinetic ,  dashed  for  hydrodynamic,) 
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ACOUSTIC  KESNaS  AS  A  rJNCTION  OF 
DmUL  MEAN  RADIUS 


Figure  11:  Kernels  as  a  Function  of  Initial  Mean  Radius.  (Solid  curve 
is  for  orthokinetic ,  dashed  for  hydrodynamic.) 
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ACOUSTIC  KERNELS  AS  A  rjNCTION  OE 
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Fi^nre  12:  Kernels  as  a  Function  of  Initial  Spectrum  Variance.  (Sol: 

curve  is  for  orth.okinetio ,  dashed  for  hydrodynamic.) 
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Figure  13:  Kernels  as  a  Function  of  Liquid  Water  Content.  (Solid  curve 
is  for  orthokinetic,  dashed  for  hydrodynamic.) 


4.2  RESULTS 


Methods  for  evsluating  the  SCE  are  discnssed  in  Prnppccher  and 

Klett  (1980)  including  the  nethod  of  Berry  and  Reinhardt  (  Berry  (1967), 

Berry  and  Reinhardt  (1974),  and  Leighton  and  Rogers  (1974)).  The  range 

-4 

of  sizes  froa  clond  droplet  to  precipitation  is  apprciimately  10  ca  to 
10  ^  cm.  In  order  to  keep  the  nnaber  of  size  classes  of  the  discrete 
spectrna  to  a  reasonable  nnaber  Berry  devised  an  exponential  subdivision 
so  that  the  drop  radius  is  given  by 

a(j)  =  a^  exp[(j-l)/J]  .  (4.6) 

In  the  present  work  the  model  code  used  is  that  of  Leighton  and  Rogers 

in  which  J  =  12/ In  2  .  There  are  130  size  classes  ranging  from 
~4  -1 

3.94x10  ca  to  6.79x10  ca.  For  each  simulated  cloud  there  is  an 

initial  droplet  density  function  as  a  function  of  droplet  mass  which  is 

specified  by  a  mean  radius  and  variance  in  a  gamma  or  Pearson  Type  III 

distribution.  The  values  selected  for  the  atmospheric  vari.;b’os  for  the 

tests  described  in  this  work  were  p  =  765  mb,  T  =  283  K,  ^  =  1. 0-2.0 
-3 

g  m 

The  model  code  was  modified  to  include  the  two  collection  kernels 
defined  by  equations  (4.3)  and  (4.4).  At  each  time  step  the  collection 
kernel  used  in  integrating  the  SCE  is  the  suas  of  equations  (4.2)  through 
(4.4) .  During  the  time  that  a  simulated  sound  is  being  applied  to  the 
cloud  the  time  step  is  one  second.  After  the  sound  is  turned  off  the 
acoustic  collection  kernels  are  set  equal  to  zero  and  the  time  step 
remains  one  second.  For  the  rest  of  the  run  the  only  process  being 
modeled  is  coalescence  due  to  gravitational  effects.  The  tests 
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discussed  Lere  are  limited  to  14  to  15  minutes  of  simulated  time  because 


numerical  instabilities  set  in  after  that.  However,  there  is  a  high 
degree  of  confidence  in  the  results  because  the  graphical  results  using 
only  the  gravitational  collection  process  have  been  compared  to  the 
results  of  Berry  and  Reinhardt  (1974b)  and  were  found  to  be  in  excellent 
agreement.  Also,  the  process  modeled  is  one  of  constant  mass  and  mass 
checks  show  gains  or  losses  of  less  than  two  percent  throughout  the 
runs . 

The  remainder  of  this  chapter  will  be  devoted  to  discussion  of  the 
model  runs.  All  of  the  results  are  presented  graphically  in  plots  that 
show  the  mass  density  functions  of  the  droplet  spectra.  The  advantage 
of  this  method  is  that  it  shows  how  the  water  mass  is  distributed  over 
droplets  of  various  radii  and  the  area  under  the  curve  between  two  radii 
represents  the  water  mass  contained  in  drops  of  those  diameters.  Also, 
in  this  type  of  ulot  equal  areas  represent  equal  amounts  of  water  mass. 
It  should  be  noted  chat  the  mean  radius  of  the  spectrum  is  the  radius  of 
the  droplet  with  mass  equal  to  the  liquid  water  content  divided  by  the 
droplet  number  concentration.  The  predominant  radius  is  that  of  the 
droplet  category  containing  the  most  mass.  The  predominant  radii  of  all 
the  cases  discussed  occur  at  the  peaks  of  the  spectra.  The 
characteristics  of  the  cases  discussed  in  the  following  pages  are 
outlined  in  Table  2. 

The  results  for  cases  1  and  2  are  shown  in  Fig.  14  through  Fig.  23. 

"3 

The  droplet  concentration  for  the  original  spectrum  is  166  cm  .  This 
spectrum  exhibits  little  change  due  to  the  effects  of  gravity  and  by  15 
minutes  shows  no  indication  of  the  development  of  precipitation  sized 
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TABLE  2 


Case 


Cloud  Spectra  Cases 


SpectruB 

Mesa  Var.  Freq. 

(it)  (Hz) 


LWC  »  1.0  g  a 


-3 


Sonnd 

lateaaity  Time 
(db)  (s) 


3 

10 

0.25 

300 

140 

5 

2 

10 

0.25 

300 

140 

10 

3 

10 

1.00 

300 

140 

5 

4 

10 

1.00 

300 

140 

10 

5 

10 

1.00 

300 

150 

5 

LWC  * 

2.0  g 

6 

10 

1.00 

300 

140 

5 

drops.  In  case  1  the  application  of  sonnd  ot  140  db  and  300  Hz  for  5 
seconds  at  the  start  of  the  ran  causes  a  shift  in  the  predooinant  radius 
and  an  increase  in  the  mass  dispersion  of  the  spectrum.  However,  even 
with  some  larger  sized  droplets  and  greater  dispersion  the  new  spectrum 
does  not  produce  precipitable  drops  in  15  minutes.  Case  2  is  similar  to 
case  1  except  the  soxind  is  applied  for  10  seconds.  The  result  is  an 
even  larger  predominant  radius  but  the  acoustically  altered  spectrum 
still  does  not  produce  precipitation  sized  drops  within  15  minutes. 

The  results  for  case  3  are  shown  in  Fig.  24  through  Fig.  28.  This 

spectrum  is  initially  more  disperse  than  in  cases  1  and  2  but  the  sound 

frequency  and  intenrity  are  the  same.  The  droplet  concentration  for 
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this  case  and  cases  4  and  5  is  225  cm  ,  Again,  the  effect  of  gravity 
alone  is  not  sufficient  to  develop  large  drops  during  the  fifteen  minute 
run.  Vhen  sonnd  is  applied  there  is  a  shift  of  mass  up  the  scale  to 
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larger  radii.  In  Fig.  28  it  is  evident  that  some  larger  drops  have 
developed  but  a  secondary  marimxim  has  yet  to  appear.  It  is  not  clear 
from  this  case  if  precipitable  drops  would  develop  if  more  time  was 
allowed  for  continued  gravitational  sedimentation. 

In  case  4  the  conditions  are  similar  to  case  3  except  the  sound  is 
applied  for  10  seconds.  At  the  end  of  the  sound  application  the  new 
spectrum  has  a  larger  predominant  radius  than  was  seen  in  case  3.  Now 
at  the  end  of  the  15  minute  run  it  can  be  seen  that  a  definite  tail  has 
formed  for  droplets  of  precipitation  size.  Berry  and  Reinhardt  (1974b) 
showed  in  their  numerical  cases  that  once  a  large  drop  signature  began 
to  develop,  mass  rapidly  moved  up  the  scale  to  the  precipitation  regime. 

The  results  for  case  5  are  shown  in  Fig. 34  through  Fig. 38.  Here 
the  effects  of  gravity  are  the  same  as  in  case  3.  The  effect  of  the 
higher  intensity  sound  is  to  shift  mass  even  further  up  the  scale  with 
an  even  smaller  mass  variance  after  the  sound  is  turned  off.  However, 
it  becomes  clear  by  ten  minutes  that  a  secondary  maximum  is  forming  for 
drops  of  precipitable  size.  The  process  is  well  under  way  by  fifteen 
minutes . 

Case  6  is  illustrated  in  Fig.  39  through  Fig.  43.  Here  the  liquid 
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water  content  has  been  increased  to  2.0  g  m  and  the  sound  parameters 
set  to  140  db  and  300  Hz  applied  for  5  seconds.  The  droplet 

_3 

concentration  is  450  cm  which  is  around  the  lower  limit  for 
continental  cumulus  clouds.  The  beginning  of  a  precipitation  tail  is 
seen  after  5  minutes  and  by  14  minutes  a  well  formed  precipitation  hump 
is  evident.  It  can  also  be  seen  that  a  precipitation  tail  begins  to 
form  by  14  minutes  due  to  gravity  alone.  However,  the  acoustic 
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agglomeration  has  produced  significant  enhancement  of  the  precipitatius 
process . 

The  model  has  been  run  vith  some  other  initial  spectra  which  will 
be  discussed  only  qualitatively  here.  When  the  initial  spectrum  has  a 
larges  mean  radius  with  the  same  liquid  water  content,  precipitation 
sired  drops  form  in  about  the  sar  i.iffle  with  or  without  the  simulated 
application  of  sound  energy.  In  cases  where  the  initial  mean  radius  is 
smaller  than  10  p  the  cur-^  at  model  does  not  initialize  correctly  so  it 
has  not  yet  been  possible  to  study  some  potentially  interesting  cloud 
spectra.  The  important  point  is  that  the  results  shown  in  this  section 
strongly  suggest  that  modification  may  take  place  for  spectra  with 
concentrations  and  liquid  water  contents  less  than  those  found  in  clouds 
over  the  continents.  Since  theory  and  experiment  show  that  acoustic 
agglomeration  is  strongly  concentration  dependent,  then  the  process  may 
be  expected  to  be  even  more  efficient  for  continental  type  spectra. 


SPECTRUM  after  5  SECONDS 


Figure  15:  Case  1  -  Spectrum  Afte',  5  Seconds.  (Solid  curve  for 

gravitational  effect  alone,  dashed  curve  for  acoustic  and 
gravitational  effects.) 
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SPECTRUM  jiJTER  5  MTSUTES 


figure  16:  Case  1  -  Spectrnm  After  5  Miimtes.  (Solid  curve  for 

gravitational  effect  alone,  dashed  curve  for  acoustic  and 
gravitational  effects.) 
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Figure  17:  Case  1  -  Spe 'trum  After  10  Minutes.  (Solid  curve  for 

gravitational  effect  aloncj  dashed  curve  for  acoustic  and 
gravitational  effect.s.) 
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Figure  18:  Case  1  -  Spectrum  After  15  Minutes.  (Solid  curve  for 

gravitational  effect  alone,  dashed  cur.e  for  acoustic  and 
gravitational  effects.) 
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Figure  20:  Case  2  -  Spectrum  After  10  Seconds.  (Solid  curve  for 

gravitational  effect  alone,  dashed  curve  for  acoustic  and 
gravitational  effects.) 
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Case  2  -  Spectrum  After  5  Minutes.  (Solid 
gravitational  effect  alone,  dashed  curve 
eravitational  effects.) 
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Figure  22:  Case  2  -  Spectrum  After  10  Minutes.  (Solid  curve  for 

gravitational  effect  alone,  dashed  curve  for  acoustic  Pud 
gravitational  effects.) 
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Figure  23:  Case  2  -  Spectrum  After  15  Minutes.  (Solid  curve  for 

gravitational  effict  alone,  dashed  curve  for  acoustic  and 
gravitational  effects.) 
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Figure  25:  Case  3  -  Spectrum  After  5  Seconds.  (Solid  curve  for 

gravitational  effect  alone,  dashed  curve  for  acoustic  and 
gravitational  effects.) 
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SPSCTRUM  AFTER  10  MINUTES 
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Figure  27:  Case  3  -  Spectrum  After  10  Minutes.  (Solid  curve  for 

gravitational  effect  alcne,  dashed  curve  for  acoustic  and 
gravitational  effects.) 
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Figure  30:  Case  4  -  Spectrum  After  10  Seconds.  (Solid  curve  for 

gravitational  effect  alone,  dashed  curve  for  acoustic  and 
gravitational  effects.) 
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Figure  31:  Case  -  Spectrnas  After  5  Minute?,.  (Solid  curve  for 

gravitf tional  effect  alone,  dashed  curve  for  acoustic  and 
gravitational  effects.) 
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Figure  33: 


''a.se  4  ~  SpectruB  After  15  Minutes.  (Solid  curve  for 
gravitational  effect  alone,  dashed  curve  for  acoustic  and 
gravitational  effects.) 
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Figure  34:  Case  5  -  Initial  Sp 
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SPECTRUM  ATtEB  5  SECONDS 
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Figure  35:  Case  5  -  Spectruia  After  5  Seconds.  (Solid  curve  for 

gravitationrl  effect  alone,  dasiied  curve  for  acoustic  and 
gravitational  effects.) 
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Figure  36:  Case  5  -  Spectrum  After  5  Minutes.  (Solid  curve  for 

gravitational  effect  alone,  dashed  curve  for  acoustic  and 
gravitational  effects.) 


SPECTRUM  JUTEE  ID  MINUTi:5 


I 

1 


R  (MICRONS) 

flO  <  10  MICRONS  150  OS 

VRR  -  I.OO  lMC  »  i.C  G/M»«3  300  ri* 


Figure  37:  Case  5  -  Spectrum  After  10  Minutes.  (Soli(3  curve  for 

gravitational  effect  alone,  dashed  curve  for  acoustic  and 
gravitational  effects.) 
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Figure  38:  Case  5  -  Spectrum  After  15  Miautes.  (Solid  curve  for 

gravitational  effect  alone,  dashed  curve  for  acoustic  and 
gravitational  effects.) 
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Figure  39:  Case  6  -  Initial  Spectrum, 
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SPECTRUM  AFTER  5  SECONDS 
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Figure  40:  Case  6  -  Spectrum  After  5  Seconds.  (Solid  curve  for 

gravitational  effect  alone,  dashed  curve  for  acoustic 
gravitational  effects.) 
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Figure  43:  Case  6  -  Spectrum  After  14  Minutes.  (Solid  curve  for 

gravitational  effect  alone,  dash'jd  curc'e  for  acorstic  and 
gravitational  effects.) 


4,3  S0U5CES  W  ERRORS 


There  are  aainly  three  potential  sources  of  error  in  the  SCE  model 
discussed  in  this  chapter.  The  first  involves  the  extrapolation  of  the 
two  acoustical  collection  kernels  to  the  conditions  associated  with 
atmospheric  clouds.  As  stated  earlier,  Chou's  numerical  model  of  the 
agglomeration  process  was  verified  by  experimental  work  but  for 
conditions  significantly  different  from  those  under  consideration  here. 
In  this  research  effort  no  experimental  data  is  yet  available  for 
verification.  Therefore,  the  formulation  of  the  SCE  model  used  here  is 
based  on  the  assumption  that  extrapolation  of  the  collection  kernels  to 
conditions  similar  to  those  in  the  atmosphere  is  valid. 

The  second  potential  source  of  error  involves  the  calculation  of 
the  collision  efficiencies  used  in  each  of  the  collection  kernels,  Tha 
use  of  the  Langmuir  formulation  in  the  orthokinetic  collection  kernel 
leads  to  collision  efficiencies  on  the  order  of  unity  for  all  the 
droplet  pairs  considered  in  the  model.  In  the  hydrodynamic  kernel  the 
use  of  the  Scott  and  Chen  formulation  leads  to  cailision  efficiencies 
which  range  from  about  .2  to  .7  for  the  variou;-.  droplet  pairs.  In 
addition,  the  coaiesence  efficiency  is  assumed  tc  equal  unity  in  all 
cases . 

Thirdly,  it  has  been  assumed  that  the  collection  volume  of  each 
collector  continues  to  be  refilled  by  some  small  droplets.  If,  in  fact, 
some  process,  such  as  condensation  or  aerosol  arift,  is  not  providing  a 
continuous  supply  of  small  droplets,  the  orthokinetic  mechanism  runs  out 
of  collectible  droplets  at  the  end  of  the  first  cycle  of  the  sound  wave. 
In  that  case  the  model  would  overestimate  the  orthokinetic  contribution 
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the  collection  process.  The  accuracy  of  these  fonanlations  will 
remain  an  unanswered  question  until  the  appropriate  erperiaentEl 
«ater  droplets  in  acoustic  fields  is  performed. 
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Chapter  V 

CONCLUDING  BEMA2ZS  AND  SUGGESTIONS  FOB  FUETHER  RESEARCH 


In  this  discassion  and  in  many  of  the  references,  theoretical  and 
experimental  work  has  been  presented  vhich  demonstrates  that  suspended 
aerosols  undertake  motion  due  to  the  propagation  of  sound  waves  through 
the  medium.  More  importantly  it  has  been  demonstrated  that,  as  a  result 
of  the  acoustically  induced  motion,  aerosols  can  be  made  to  agglomerate. 
Under  the  right  conditions  an  aerosol  spectrum  can  be  significantly 
modified  and  precipitation  of  ths  aerosol  initiated.  This  agglomeration 
process  can  be  very  efficient  for  the  concentrations  and  dispersions 
associated  with  industrial  pollutants.  The  purpose  of  the  work 
described  in  this  thesis  has  been  an  initial  examination  of  the 
possibility  that  acoustic  agglomeration  can  also  be  an  efficient  process 
in  the  atmosphere.  A  final  determination  on  the  matter  will  not  be  made 
until  the  experiments  are  ac.ually  performed  on  clouds.  However,  the 
results  of  Chapter  IV  indicate  that  if  the  physical  processes  leading  to 
the  agglomeration  of  water  droplets  in  a  sound  field  are  adequately 
represented  in  the  numerical  model  used  here,  then  acoustic 
agglomeration  can  be  an  efficient  process  in  the  atmosphere. 

The  first  extension  of  this  work  should  be  to  modify  the  SCE  model 
to  eliminate  the  numerical  instabilities  that  limit  the  liquid  water 
contents  anJ  lengths  of  the  model  runs.  This  can  be  accomplished  by 
revising  the  code  to  the  original  specifications  of  the  Berry  and 
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Eoinhardt  model.  The  second  extension  should  be  to  modify  the  code  of 
the  model  to  alloi?  computation  of  the  initisi  spectrum  fsos- 
specif ication  of  the  cloud  liquid  water  content  and  droplet 
coaeeatrstion.  To  numeiicslly  examine  the  effect  of  acoustic 
agglomeration  on  continental  cumulus  clouds  will  require  iatial iration 
with  spectra  that  are  narrower  with  higher  concentrations.  These  clouds 

-3 

typically  ha've  concentrations  of  500-1500  cm  with  liquid  water 

-3 

contents  of  1. 0-3.0  g  m 

There  are  two  areas  of  meteorology  in  which  acoustic  agglomeration 
may  play  a  significant  role.  The  first  is  as  a  mechanism  for  artificial 
weather  modification.  Experimental  work  should  be  designed  to  test 
acoustic  agglomeration  of  water  droplets  in  a  cloud  chamber,  if 

possible,  or  under  field  conditions.  This  would  help  verify  the 

extrapolation  of  the  acoustic  collection  kernels  used  in  this  modeling 
study.  It  would  also  provide  insight  into  the  possibly  adverse  effects 
of  acoustically  induced  evaporation,  an  effect  which  has  been  reglected 
in  this  study.  Additionally,  field  studies  could  provide  information  on 
the  distance  limitations  imposed  by  the  attenuation  of  sound  in  air 
wheieby  the  high  intensities  required  for  agglomeration  are  reducid. 

The  second  area  of  interest  is  the  study  of  the  effect  of  thunder 
on  clouds.  In  addition  to  studies  of  the  fundamentals  of  cloud  droplet 
interactions  with  acoustic  energy  and  the  exploration  of  the  technical 
fea«ibility  of  artificial  applications,  the  role  of  thunder  as  an 

augmentor  of  the  natural  precipitation  process  deserves  attention. 

Extension  of  the  piesent  work  to  the  study  of  thunaer  should  include 
expanding  the  model  to  simulate  a  multi-frequency  sound  spectrum  similar 
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to  that  produced  by  thunder.  Also,  since  thunder  is  initiated  as 
shoch  wave,  investigation  of  that  type  of  wave 


behavior 


should  b 
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